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Red giant branch stars as probes of stellar populations. II. 
Properties of the newly discovered globular cluster GLIMPSE-C01 
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Abstract. Deep near infrared photometry of the newly discovered Galactic globular Cluster GLIMPSE-C01 is reported. We 
derived for the first time the metal abundance of this object from the slope of the RGB: [Fe/H]=- 1.61 ±0.14 in the scale of Zinn 
(as implemented in Harris [T99oV [Fe/H]=-1.44±0.12 in the scale of Caretta & Gratton 119971 . and [Fe/H]=-1.12±0.12 in the 
scale of Ferraro et al. ( 1999 1. The tip and the clump of the red giant branch were used to confirm the estimates of Kobulnicky 
et al. y005il, placing the cluster at D~3.7±0.8 kpc, behind Ai/~15 mag of visual extinction. The best fit to the radial surface 
brightness profile with a single-mass King's model 1 1962 1 yielded core radius r c =0.78 arcmin, tidal radius r,=27 arcmin, and 
central concentration c=1.54. 

Finally, we estimate the number of the "missing" globulars in the central region of the Milky Way. Based on the spatial dis- 
tribution of the known clusters, and assuming radial symmetry around the Galactic center, we conclude that the Milky Way 
contains at least 10±3 undiscovered objects. The distribution of known clusters in the bulge seem to resemble the orientation 
of the Milky Way bar. 
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1. Introduction 

The all-sky infrared surveys carried out during the recent years 
have brought the discovery of a number of new clusters, hid- 
den by the dust extinction in the plane of the Milky Way. They 
usually suffer Ay>10-20 mag of extinction, making them in- 
visible in the optical wavebands. The vast majority of these ob- 
jects appear to be a few million years old (Ivanov et al. 2002 
2005] Borissova et al. l2"0"0"3l I20TJ5I but a few have proved to 
be analogues of "classical" globular clusters (Hurt et al. 2000 
Kobulnicky et al. EOOSl Carraro l2"0Tj5V 

The globular cluster GLIMPSE-CO 1 was discovered 
(Kobulnicky et al. 2005 1 during the Galactic Legacy Infrared 
Mid-Plane Survey Extraordinaire (hereafter GLIMPSE; 
Benjamin et al. 2003 1). The survey is mapping the Galactic 
plane from |1| = 10° to 65° and |b|<l° with the Infrared-Array 
Camera (IRAC; Fazio et al. 2004) on Spitzer Space Telescope 
in 3.6, 4.5, 5.8, and 8.0 micron bands. It was also identified in- 



dependently by Simpson & Cotera ( 2004 1 who cross-correlated 
ASCA X-ray and IRAS sources with the 2 MASS. 

The isochrone comparison indicated that GLIMPSE-CO 1 
is indeed a globular cluster, at least a few Gyr old. The near 
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infrared color-magnitude diagram (CMD) suggested that the 
object suffers Ay=15+3 mag of extinction, and the analysis of 
the line-of-sight 13 CO yielded a distance of 3.1-5.2 Kpc, based 
on a cinematic model of the Milky Way. The question of the 
cluster abundance remained open, as well as the need to verify 
the distance. 

Here we report properties - metallicity, extinction, and 
distance - derived from deep near infrared photometry of 
GLIMPSE-CO 1. We also estimate the possible number of 
undiscovered globular clusters in the central region of the 
Milky Way. 

2. Observations and data reduction 

The JHKs imaging observations of GLIMPSE-CO 1 were car- 
ried out in Nov 2004 under non-photometric condition and 1 
arcsec seeing with the Sofl (Son of ISAAC) at the NTT. The 
instrument is equipped with Hawaii HgCdTe 1024x1024 de- 
tector, with pixel scale of 0.288 arcsec px _1 . 

We took 12 images in each filter, in jittering mode with 3 
arcmin jitter box size to ensure that on each image the cluster 
is located on a different place on the array. Each image was 
the average of 3 frames of 10 sec for J, 13 frames of 3 sec 
for H, and 20 frames of 3 sec for Ks , comprising total integra- 
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Fig. 1. True-color near infrared composite image of GLIMPSE- 
CO 1. The J, H, and Kg bands correspond to blue, green and 
red, respectively. North is up, and East is to the left. The field 
of view is 4.92 arcmin on the side. The dark lane in the North- 
West part of the image corresponds to a region of mid-infrared 
emission, and it is a real feature. (This figure is available in 
color from the electronic edition of the journal.) 

tion times of 6, 7.8, and 12 min, respectively. Individual images 
were sky-subtracted, flat-fielded, aligned, and combined into a 
single image. 

The stellar photometry of the final images was carried out 
with ALLSTAR in DAOPHOT II (Stetson [T9931 . The typical 
photometric errors vary from 0.03 mag for the Ks~10 mag 
stars to 0.10 mag for Ks~13 and 0.15 mag for Ks~16 mag. 
The photometry of the the faint stars is affected by confusion, 
especially near the cluster center but our metallicity estimates 
is based on stars with Ks<13 mag. 

The photometric calibration was performed by comparing 
our instrumental magnitudes with the 2 MASS measurements 
of 67 stars, covering the color ranges 0.61<./— flf,s<5.98 and 
0.21<H-K S <2.32 mag, and magnitude range 8.33<.K S <13.13 
mag. The transformation equations are: 

J-K s = (0.981 + 0.005) x (j -k) + (3.208 + 0.007), (1) 

H-K s = (0.971 + 0.016) x (h - k) + (1.820 + 0.014), (2) 

and 

K s = (1.008+0.007)xfc+(0.0126±0.006)x(7-^ s )-(10.096±0. 

with r.m.s. =0.047, 0.062, and 0.044 mag, respectively. Here k, 
j-k, h-k are the instrumental magnitudes and colors and Ks , 
J-K$, H-Ks are the magnitudes and colors in the 2 MASS 
system. The standard error values are given after each coeffi- 
cient. 

A true-color composite of GLIMPSE-CO 1 is shown in 
Figure H The image indicates substantial variations of the ex- 
tinction in the field. Dust structures are present even across the 
cluster face, suggesting that the completeness of the photom- 
etry varies significantly. The crowding also contributes to this 
effect, as it can be seen from the luminosity function of stars in 
the field of GLIMPSE-CO 1 (FigureEJ: it is clear that the com- 
pleteness of the photometry in the cluster area (long-dashed 
line) is at least a magnitude shallowed than that in the field 
(dotted line). 

3. Properties of GLIMPSE-C01 

3.1. Structural parameters 

The structural parameters of GLIMPSE-CO 1 were determined 
using the iterative star counts method of King ( 1962 1 after ran- 
domly removing the field stars (see Sec. l3.3> . The radial profile 
was built only from the stars in the South-Eastern half of the 
GLIMPSE-CO 1 to exclude the effects of the variable extinction 
across the cluster face. The best fit to the radial surface bright- 
ness profile with a single-mass King's model (FigurefJ} yielded 
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Fig. 2. Near infrared luminosity function for: all stars in the 
field of GLIMPSE-CO 1 (solid line), stars within R<30 arcsec 
from the center of the cluster (long-dashed line), and stars in an 
annulus with inner radius of 50 arcsec and the same area as the 
cluster region (dotted line). 



.o^re^adius r c =0.78 arcmin, tidal radius r f =27 arcmin, and cen- 
tral concentration c=1.54. The inhomogeneous foreground ex- 
tinction - significant even in the near-infrared - makes it im- 
possible to confirm the suggestion of Kobulnicky et al. (|2005 1 
that GLIMPSE-CO 1 is elongated but we certainly can not ex- 
clude such a possibility. 

3.2. Distance and extinction 

The CMD of GLIMPSE-CO 1 is plotted in Figure H The left 
panel contains all stars with J and Ks photometry, and shows 
the presence of a few sequences. The nature of these sequences 
become apparent from the other two panels: the cluster (middle 
panel) shows a red giant branch (hereafter RGB) at J-Ks ~3-4 
mag, and the field (right panel) is dominated by main sequence 
at J—K s ~\-2 mag and a group of red stars with J-K s >3 mag, 
possible just highly reddened background objects. 

An inspection of the CMD reveals the presence of a clump 
of stars in the red giant at A^~13.0 mag. It is also noticeable 
on the cluster luminosity function (Figure [2j. Taking into ac- 
count the width of the structure along the K s axis, we assign an 
uncertainty of 0.25 mag (equal to the bin size used to build the 
luminosity function) to the apparent magnitude of the clump. If 
this structure is indeed the red clump, it has an absolute mag- 
nitude of M/j-=-1.61+0.03 mag (Alves 2000l, yielding a dis- 
tance modulus {m-M)k~ 14.6+0.3 mag. Note that this calibra- 
tion doesn't take into account any metallicity effects (see for 
discussion Salaris & Girardi 2002 1. 
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Fig. 4. Near infrared color-magnitude diagram for: all stars in the field of GLIMPSE-CO 1 (left panel), and for starts in two regions 
with identical areas - a circle with radius R=20 arcsec around the cluster center (middle panel), and an annulus around the cluster 
with inner radius of 60 arcsec (right panel). The numbers of stars plotted in each panel are given in brackets. In the middle panel: 
the best fit to the RGB is drawn with a solid line; crosses mark stars included in one realization (see Section l3~3l for explanation). 



The distance to the cluster can also be measured using the 
RGB tip (i.e. Ivanov & Borissova 2002 1 but this method suf- 
fers from strong metallicity dependence and it is hampered 
by the small number of stars at the upper end of the RGB. 
Nevertheless, we carried out this test as a consistency check. 
If the tip is located at Ks~8.7 mag, and assuming Mk=-6.1 
mag for [Fe/H]=-1.5 (Ivanov & Borissova 2002 1 we obtain 
(m-M)^~14.8 mag, in agreement with the estimate given 
above. 



The color-color diagram of the cluster field is shown in 
Figure|5] Our data confirm the visual extinction Ay~15+3 mag 
derived by Kobulnicky et al. (2005 1. Throughout this paper we 
used the reddening law of Rieke and Lebofsky ( 1985 1, giving 
Ajf~ 1.7+0.3 mag. This yields reddening-free distance modu- 
lus of (m-M)o~\2. 9+0.4 mag, or D~3. 8+0.7 kpc, in agree- 
ment with Kobulnicky et al. J2005> . The CMD shows a no- 
table spread of colors among the fainter stars that might be 
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Fig. 3. Structural parameters of the cluster: King's profile fit. 

contributed to differential extinction and to contamination from 
extended emission and background sources. 

3.3. Metal abundance 

The RGB slope allows to derive abundances of globular 
clusters (i.e. De Costa & Armandroff 1990 1 because of the 
metallicity-dependent opacities in cool stars. The RGB slope is 
independent of the reddening, which is an important advantage 
in studies of heavily obscured clusters such as GLIMPSE-CO 1. 
Here we applied the analysis and the calibrations, described in 
Ivanov & Borissova J2002> . 

The first step was to remove the fore- and background stars, 
which constituted ~8-9% of the stars in the designated cluster 
area - a circle with 40 arcsec diameter. We used a Monte-Carlo 
technique: the CMD of the cluster was divided into rectangu- 
lar bins. Similarly, the CMD of a circular annulus, centered on 
the cluster, with an inner radius 60 arcsec, and the same are 
as the cluster region. Then, from each bin of the cluster CMD 
we removed randomly a number of star, equal to the number of 
stars in the corresponding bin of the "field" CMD, producing a 
single realization. Next, we carried out a least square fit to the 
RGB (Ivanov & Borissova 2002 see also Appendix lAl in this 
work): 



(J-K s ) =ZP + SlopexK s 



(4) 



on the corrected diagram to derive the slope and the zero 
point. The RG locus was defined after inspecting the CMD: 
8.5</^s<13.0 mag, and 2.5<J-Ks<4.0 mag. The bin-size was 
0.5 mag along both axis. We removed the 10-<x outliers, and 
repeated the fitting. Typically, the fitting coefficients obtained 




Fig. 5. Near infrared color-color diagram for starts in two re- 
gions in the vicinity of GLIMPSE-CO 1, with identical areas - a 
circle with radius R=30 arcsec around the cluster center (open 
circles; the size is proportional to the apparent brightness of 
the star), and an annulus around the cluster with inner radius if 
50 arcsec (solid dots). A vector representing visual extinction 
Ay=15 mag is plotted. The colors of RGB (solid line) and the 
main sequence stars (dashed line) are also shown (Frogel et al. 
1978 no correction for the photometric systems is applied). 



in the two iterations were statistically indistinguishable. Note 
that the spread of colors along the RGB is dominated by differ- 
ential extinction, rather than measurement errors. Total of 10 4 
realizations were obtained, and we averaged the fitting results. 
The distributions of the RGB slope and the zero point of the fits 
are shown in Figure [6] 

Finally, we derived the metal abundance from the 
calibration of Ivanov & Borissova J2002> . and obtained 
[Fe/H]=-1.61+0.14 in the scale of Zinn (as implemented in 
Harris [1996), [Fe/H]=- 1.44+0. 12 in the scale of Caretta & 
Gratton iTW7\ . and [Fe/H]=-1.12+0.12in the scale of Ferraro 
et al. J1999I . The calibrations of Valenti, Ferraro, & Origlia 



(2004) wield [Fe/H]=- 1.14+0. 16 and [Fe/H]=-0.97±0.15, for 
the last two scales. Both values are systematically lower than 
ours but the differences are within 1-2 <x. 

These measurements rely on the important assumption that 
GLIMPSE-CO 1 has the same age as the other Milky Way glob- 
ular clusters, used to derive the metallicity versus the RGB 
slope calibrations. The comparison of the cluster CMD with 



stellar isochrones (Kobulnicky et al. 2005 ) seems to indicate 
an age of at least 8 Gyr but this issue can not be addressed until 
more accurate age estimate becomes available. 
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Fig. 6. Distributions of the RGB slope (top panel) and the zero 
point (bottom panel) of the fits for the 10 4 Monte-Carlo re- 
alizations with random removal of the fore- and background 
contamination. Gaussian fits to the distributions are overplot- 
ted, and the average values with their standard errors are given. 
See Section EOl for details. 



4. The case of the missing clusters 

The discoveries of GLIMPSE-CO 1 and Whitting- 1 prompted us 
to reinvestigate the question if there are any more undiscovered 
globular clusters in the Milky Way. The total number of Milky 
Way globulars estimated by van den Bergh (1998 1 is 160+20, 
slightly above the currently known ~ 150 objects. Barbuy, Bica, 
& Ortolani ( 1998, see their Figure 3 and 4) argued that some 
missing clusters are probably located in the general direction 
of the Galactic Center. Here we adopt a modification of their 
method. The XYZ Galactic coordinates used in this section are 
taken from Feb 2003 edition of the Milky Way Globular Cluster 
Catalog (Harris 1996 1 with their face value. A possible caveat 
is the variation of the ratio of the total extinction to the redden- 
ing Ry. For example, a change or Ry from 3 to 3.3 introduces 
10% change in the distances, or almost 1 kpc at the distance 
of the Galactic Center. The averaging over different dust prop- 
erties along the line of sight helps to minimize the related dis- 
tance uncertainties but the variations of Ry toward the bulge 
definitely need further investigation. 

First, we assumed that all missing clusters are located in 
the Galactic plane and close to the Galactic Center. Therefore, 
we consider only the region with |Z|<0.5 kpc and Rgc^3.0 kpc. 
These constraints reflect the spatial distribution of the obscur- 
ing material. They make our estimate of the missing clusters 
only a lower limit, as it was demonstrated by the recent discov- 
ery of the off-the-plane cluster Whiting- 1. 

An observer at the Galactic Center should detect an equal 
number of globulars toward the Sun, in the antisolar direction, 
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Fig. 7. Spatial location of the globular clusters in the selected 
region (see Section |4} projected on the Galactic plane. The 
Sun is marked with S, the Galactic Center lies at the origin 
of the plot. The X-axis increases toward the Sun, and the Y- 
axis increases in the direction of Galactic rotation. The lo- 
cations of 2MASSGC01, 2MASSGC02, and GLIMPSEC01 
are indicated (2M1, 2M2, and GL1, respectively). The dot- 
ted line shows the position of the Galactic bar, according to 
Weinberg ( 1992 1. Dashed lines indicate the limits of the bins 
for the histogram shown in the inset - it gives the distribution 
of the number of clusters in the selected region viewed from 
the Galactic Center. The viewing angle is deg in the direction 
toward the Sun and it increases to 180 deg in the antisolar di- 
rection. Poisson uncertainties are shown with bars. The area of 
the shaded zone limited by the average level of the outer and 
the inner bins corresponds to the number of the clusters, needed 
to obtain a homogeneous distribution. 

and in the directions perpendicular to the line, connecting the 
Sun and the Galactic Center, i.e. the globular cluster distribu- 
tion along any direction in the Galactic plane should be flat. 
However, it show some structure (Figure [7J: there are more 
clusters in the directions toward the Sun and in the opposite 
direction than in the direction perpendicular to the line con- 
nected the Galactic Center and the Sun. Ten additional clusters 
are necessary to flatten the histogram (indicated by the shaded 
area), setting a lower limit of the missing clusters. 

This estimate has to be treated with caution for a number 
of reasons. For example, it is sensitive to the adopted globu- 
lar cluster distance, and may contain systematic errors related 
to the adopted globular cluster scale (and the Galactic Center 
distance). In addition, while the location of the clusters on the 
sky are known very well, their distances contain statistical ob- 
servational errors, artificially extending the distribution along 
the line of sight, similarly to the "finger of God" effect known 
to the extragalactic astronomers. This would move some of the 
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clusters from the middle two bins to the outer bins, increasing 
estimate. 

Finally, the spatial location of the clusters may be affected 
by the presence of a bar or a triaxial bulge in our Galaxy. The 
major axis of the bar has an angle of only about 36+ 10 deg from 
the line of sight toward the Galactic center (Weinberg 1992 1, 
with the near end lying in the first quadrant. The major axis of 
the bar is close to the direction in which the cluster distribution 
is elongated. 

Our result suggests that the searches for hidden clusters 
carried out so far are biased toward the regions closer to the 
Galactic Center so the further away from the Galactic Center is 
a region, the less likely it is to get attention. This implies that 
the future searches has to modify their strategy, including re- 
gions in the Galactic plane, as far as 7-15 deg from the Galactic 
Center. 

5. Summary 

We report deep near infrared photometry of the newly dis- 
covered Galactic globular Cluster GLIMPSE-CO 1, and we de- 
rived for the first time the metal abundance of this object from 
the slope of the RGB. The cluster appears metal-poor, with 
[Fe/H]=-1.6 in the scale of Zinn. We confirm the distance and 
reddening estimates of Kobulnicky et al. (2005), placing the 
cluster at D~3.7+0.8 kpc, behind Ay~15 mag of visual extinc- 
tion. 

Finally, we estimate the number of the missing clusters in 
the central region of the Milky Way. Based on the location of 
the known clusters, and assuming radial symmetry of the clus- 
ter distribution around the Galactic center, we conclude that the 
Milky Way contains at least 10+3 undiscovered globular clus- 
ters. 
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Appendix A: Red Giant Slope - Metallicity 
calibration 

The equations given bellow superseded Eqns. 1-4 from Ivanov 
& Borissova ( 2002 1 which contained an error. 

(J-K s )o=RGB ZP + RGBsixM Ks (A.l) 

RGBs, = a s {) ' + of x [Fe/H] (A.2) 

RGB ZP = al P + of x [Fe/H] (A.3) 

[Fe/H] = [(J-K s )o-a s 'xM Ks -of] / [a\ p xM Ks +af ](A.4) 
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